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Abstract 

The problem of velocity selection of reaction-diffusion fronts has been widely investigated. While 
the mean field limit results are well known theoretically, there is a lack of analytic progress in those 
cases in which fiuctuations are to be taken into account. Here, we construct an analytic theory 
connecting the first principles of the reaction-diffusion process to an effective equation of motion via 
field-theoretic arguments, and we arrive at the results already confirmed by numerical simulations. 
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1 



I. INTRODUCTION 



Reaction-diffusion front propagation in nonequilibrium systems is a topic that has been 
receiving an increasing attention recently. The numerous possible applications of the theory, 
as can be systems like flames bacterial colonies or population genetics 0], is of 
course one of the reasons of this recent interest. One of the most common approaches to 
this problem has been the use of deterministic reaction-diffusion equations, like the Fisher 
equation Q]. This equation, that combines logistic growth with diffusion, is one of the most 
important mathematical models in biology and ecology [4]. In one spatial dimension, the 
Fisher equation reads 

dtU = Dd^^U + aU - bU^. (1) 

One can think this equation as the mean field description of a reaction-diffusion process of a 
single species of random walkers A undergoing the reactions of birth A ^ A+ A at rate a and 
annihilation A+A — at rate 6/2. The analysis of this equation is straightforward. Consider 
the boundary conditions U ^ b/a when x — oo and f/ — when x oo. Thus the 
linearly stable phase b/a invades the linearly unstable phase 0. Assuming a stationary front 



profile U(x,t)=U(x-vt)=U(z) and shifting variables x y'D/ax, t t/a and U {a/b)U 
we get 

U" + cU' + U-U^ = 0, (2) 

where c = v/ \J Da. The velocity of the front is controlled by its edge, this means, the region 
of the front that is closer to the unstable phase U = 0. We can thus linearize Eq.® around 
this value to get 

U" + cU' + U = 0. (3) 
The only physically acceptable solution to this equation is 

U{z) ~ e-^^ (4) 

and substituting Eq.(jH) into Eq.Q we get that 

c = 7 + -, (5) 

7 

for an arbitrary 7. It is clear that the range of velocities is thus c > 2, and it was shown 
that the minimal velocity is selected in the long time limit P]. We can thus conclude that, 
in this limit, v = 2\/ Da. 



We will now show that this picture changes strongly when internal fluctuations effects, 
due to the flnitness and discretness of the reactants, are taken into account. 



II. THE FIELD THEORY 



We will consider a single species particles A undergoing the reactions A ^ A + A at rate 
a and A + A ^ ^ at rate A. Further, we suppose the particles A performing a random walk 
in a one dimensional lattice with lattice spacing b. The exact description of the problem is 
given by the following master equation: 
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with 
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= ^X][(^e + l)P{...,ni - l,ne + 1, - niP{...,ni,ne, ...;t)], (7) 

D {e} 



where {e} denotes the set of nearest-neighbor sites adjacent to i and D is the diffusion 
constant, 



dP{{niy,t) 



dt 



(7[{ni - l)P(...,ni - 1, - niP{...,ni, 



(8) 



and 



dP{{ni};t) 



dt 



\[{ni + 2){ni + l)Pi...,ni + 2,...;t)-niini-l)Pi...,ni,...;t)]. (9) 



For simplicity we will choose an uncorrelated Poisson distribution as initial condition for 
our master equation: 

rii 

P(W;t = 0) = e-^(°)n^, (10) 



where N{0) = X^j'^oi- We can map this master equation description of the system into 
a quantum fleld-theoretic problem. This connection was flrst proposed by Doi 6], further 
elucidated by Peliti Q] and a deep generalization of it can be found in the influencing article 
by Cardy and Tauber j^. We can write this theory in terms of the second-quantized bosonic 
operators: 



[tti, ttj] = 0, [al, a]] = 0, tti |0) = 0, 



(11) 
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whose effect is to create or to anniliilate particles at tlie corresponding lattice site: 

4 \--,ni, ...) = \...,ni + 1, ...) , (12) 

tti \...,ni, ...) = Hi \...,ni - 1, ...) , (13) 
where we have defined the states as: 

iK}) = n(«lno)- (14) 

i 

Thus we can define the time-dependent state vector as: 

{Ui} 



and claim that it obeys the imaginary time Schrodinger equation 

d 
di 

with the hamiltonian 



i\m) = -H\m), (16) 



^ = J2y-§Yl «^^(«- - - - (^«^)']^*' + - 4Wi(^ij ■ (17) 

Note that we recover Eq.® if we substitute Eq. (|TK|l and Eq. (|T7jl in Eq. (|T!)|) . The time 
dependent expectation value of an observable O is given by: 

{0{t)) = J20{{n.})P{{n,};t). (18) 

{rii} 

To compute this quantity in the field-theoretic formalism we need to introduce the Glauber 
state: 

{S\ = {0\l[e'^% {S\0) = 0. (19) 

i 

Note that this state is a left eigenstate of the creation operator with eigenvalue 1, implying 
that for any normal-ordered polynomial of the ladder operators one has 

{S\Q{{al},{a,}) = {S\Q{{l},{a,}). (20) 

Thus we can write expectation value Eq.()18|) as 

{0{t)) = {S\0{{a,})\m)- (21) 



We can write this expectation value as coherent-state path integrah 



(0(T)) 

where the action is given by 



S[i,„i^f,T] = J2(£dt 



Performing the continuum hmit: 



If 

Y,[i'e{t)-i^m^y'-^ij{x,t), 

{e} 



we get the action: 



S[^l}, il)]T]= I dx 



dt ip{x,t) 



dt dx'^ 



ij{x,t) 



-Ao[l - ipix, tY]ip{x, tf + a[l - i){x, t)]i){x, t)ip{x, t) 



where Xn = b\. 



(22) 



(23) 



(24) 



(25) 



III. PERTURBATION THEORY 

In order to study perturbatively this field theory we will perform a change of variables 
to rend the action Eq. lj^ dimensionless: 



t — >■ — , X 
a 



-X, 1p ^ Ipy 1p ^ t4' 

A 



(26) 



this way we get: 



S {ip , ip) = e / dxdt [ ilj{x,t) 



d 



dt dx"^ 



il){x,t) 



[1 - ^(x,t)2]^(x,t)2 + [1 - V'(x,t)]V'(x,t)^(x,t) 



(27) 



where e 



-1 _ -/D^ 



, -. We will use from now on some standard results involving functionals 

AO n 

and functional integrals in field theory, they can be seen for instance in The functional 
Z{fj,ri) with external sources is: 

Z{fi,ri) = /"DV'(a:,t)^V'(a;,i)e"^^^+^'^'^'^*f^^^'*^'^^'^'*^+'^^'''*^''^'^'*^'). (28) 



Using the steepest-descent procedure, we know that the functional integral Eq.(j2HI) in the 
limit e ^ is dominated by the saddle points: 

ss 

fi{x,t), (29) 
riix,t). (30) 



SS 



5%l){x, t) 

Eqs. fl29l30p in the absence of the external sources (r^ = = 0) are the mean- field equations for 
the reaction-diffusion process. We will study perturbatively the functional integral Eq. ()28|) 
in a neighborhood of the "classical field", say, the solutions of the saddle-point equations 
(I29I3U|) :'?/'^. i)... Thus we will use the expansion: 



and expanding the action in powers of e we find: 



(31) 
(32) 



Si^.ijj) - r]^! - fjip = S{ipc,ipc) - #c - fjipc + 

dXidX2dtidt2 ir-T TTT-T TT Xl^^l, il)x(a;2, ^2) + 



6ilj{xi,ti)6ilj{x2,t2) 
6^S 



ip=i'c 



5il){xi,ti)5i>{x2-,t2 



X{Xl,ti)x{x2,t2) + 



6^S 



5i}{xi,ti)5i){x2,t2 



(33) 



X{.Xi,ti)x{,X2,t2) 

It is very important to note at this point that it is the edge of the front that leads to the 
marginal stability criterium, say, to the velocity selection. And the edge of the front is 
characterized by a low occupation number, so we can neglect the terms proportional to 
and x^, that reflect the presence of more than one particle at the corresponding site, as we 
consider this event to be unlikely if we go far enough in the edge. We can see this clearly 
if we remind that x S'-nd x ^''^^ the eigenvalues of the annihilation and creation operators 
respectively, and this way any of them squared reflects the possible presence of two particles 
in the same place. 



The functional integral at this order becomes: 

Z{ri,T]) ^ Zo{r],r]) DxDxe «<^c(^i,ti)<5i/.c(^2,t2) , 



(34) 



where 



and therefore: 



-:[S{i>c,ipc)-fltpc-i'cri] 



Z{ri,ri) = NZo{ri,ri) 



det- 



6^S 



-1 



5llJc{Xi,ti)6llJc{x2,t2) 

The normahzation factor is fixed by the condition Z{0,0) = 1. 

The connected generating functional W{fj,ri) = elnZ{'fj,ri) at this order is then: 



(35) 



(36) 



W{ri, ri) = Woiv, r]) + eWiir], r]) + 0(e' 



where 



tr In 



6^S 



— tr In 



5^S 



f]=ri=0 



(37) 



(38) 



5l})c{Xi,ti)6l})c{x2,t2) 

Let us now perform the Legendre transformation 

r(0, 4>) = j dxdt[r]{x, t)4>{x, t) + r)(x, t)0(x, t) - ^0(17, r/) - eWi{f], r/)] + 0{t^). (39) 
At one-loop order the IPI functional is: 

r(0, 0) = 5(0, 0) + eri(0, 0) + 0(e2), (40) 

with: 



ri( 



tr 



In- 



5^S 



-In 



0=</.=O. 



(41) 



50(xi, ti)50(a;2, ^2) 5050 
We can interpret the IPI functional as an effective action that will lead us to new effective 

equations of motion: 



5T 



5T 



0. 



In our particular case, action (P7|) gives: 

5^S 



(^i- a,.,. + 1)^2 + (400 -20)^2, 



(42) 
(43) 

(44) 



where 82 = 5{xi — X2)5{ti — 12). We conclude therefore: 



ri( 



dxdt (x, t \ In 



l + idt- d,, + 1)^^(400 - 20)J \x, t) = 
/d.*(MlB-4. + l)-'(4*-2,)|M). 



(45) 



where we have made use of the low occupation number approximation for the front edge 
and a power series expansion of the logarithm. If we evaluate the matrix element we get: 

Ti = C y dxdt{4:4)(t) - 20), (46) 

where the constant C is given by: 

j dp j dte{t)e~^P'-^^^' j dwe'""' 

where we have used the fact that 0(0) = 1/2 (this property of the Heaviside G function 
can be found more rigourosly proven in [^) and that the integral over the whole momentum 
space is the volume of the first Brillouin zone, where ho is the dimensionless lattice spacing. 
Thus the effective action reads: 

r(0, 0) = ^(0, 0) + e-^ / dxdt{A^<P - 20) + 0(e2), (48) 
2oo J 



and the new equations of motion Eqs. ()42|4Hj) are: 

{dt - a,,,,)0 + 2002 - 200 + + e-^(40 - 2) = 0, (49) 
-{dt + 9,.)0 - 2(1 - 02)0 + (1 _ 0)0 + e-^(40) = 0. (50) 

If we remind that the field represents the expected value of the front density, and that 
we are everywhere supposing that we are on the front edge, this quantity should be small 
enough to consider e0 neligible. This way we get that = 1 + 2e/6o solves equation Eq. ljMUl . 
Substituting this result in Eq. ()49|) and taking into account that 6o = \fojDh we get 

9t0 = 5..0 + 0-202 + -. (51) 

a 

Note that in this case, contrary to the mean field approach, a positive phase propagates 
into an (infinitesimaly) negative phase, something that is clearly unphysical. This is the 
deterministic expression of the compact support property of the front, i.e., the front becomes 
identically zero at a finite value of x. Actually, this property has been rigourosly proven for 
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this kind of fronts Q|. In this regime, since there are no particles, no reaction is possible, 
only diffusion from adjacent sites is allowed, leading to the action: 



S, 



diff 



JdxJ dt[4>{dt - d^^)<j)]. (52) 



Since this action is quadratic, it is the effective action to any order, and this implies that 
the effective equation of motion for the front propagation at first order in A/cr is 

dt^ = + (^0 - + ^ 6(0). (53) 

This last derivation deserves a further explanation. It may be surprising to the reader that 
Eq. ljSlj) performs such an unphysical behaviour, but it is actually what one would expect a 
priori. Indeed, the analysis of an effective action commonly yields a shift of the fixed points 
of the original one, that is what has happened here. This suggests that the correct physical 
interpretation of the problem should have been the corresponding to a moving boundary 
one. This is, at the beginning, we should have had into account two different actions, one 
for the space full of particles and one for the empty space, and study the propagation of the 
boundary between them. This preserves the physical meaning all along the derivation. It 
might be desirable to solve this problem without splitting it into t wo p arts, something that 
maybe could be done by using stochastic differential equations 



;0 t wo P e 

HQ. 



IV. FRONT PROPAGATION AND VELOCITY SELECTION 

To study how the front propagates let us perform the change of variables = m — A/cr. 
At leading order, Eq. (j^ becomes: 



dtu = dxxU + {u — 2u^)Q — — ^ . 



(54) 



Clearly, fields u and propagate at the same speed. Eq. (jMj) was heuristically proposed and 
studied by Brunet and Derrida Q], and we will summarize the main conclusions of their 
work. We will consider that for sufficiently long times the front will converge to a stationary 
shape, this is, u{x,t) = u{x — ct) = u{z), and Eq.(jnH) becomes 

u" + cu' + {u - 2u'^)e (u - = 0, (55) 
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where c is the dimensionless front speed. We can distinguish between three regions in the 
front edge, in the first one, the equation of motion is given by 

u" + CU +u- 2u^ = 0, (56) 

in the second one, the field is small enough that we can linearize to get: 

u" + cu' + u = 0, (57) 

and in the third one u < A/a, so we get 

u" + cu = 0. (58) 

We impose as boundary conditions the continuity of the first derivative between the bound- 
aries of the three different regions. It can be shown that this leads to a dimensionless 
velocity 

at first order in X/a. The front velocity with the corresponding dimensions is thus: 



( 2 - T-27TTT ) • (60) 



This correction has already been confirmed in numerical simulations 
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la UM, and it 



shows a very slow convergence to the mean-field velocity in the limit A/cr — » 0. This shows 
that the discretness of the reaction process strongly shifts the velocity to a slower one. 

As a final remark, we would like to underline that the cutoff derived is not a particularity 
of the low dimensional topology of the problem. Indeed, if one considers the (i-dimensional 
problem and performs all the calculations shown here for the particular case d = 1, one 
arrives at the dimensionless equation: 

dt^ = V'4>+ (4>-24>' + -]q{^), (61) 



a ^ 

locally describing the edge of the front. The reason of this independence between the ap- 
pearance of the cutoff in the reaction term and the dimensionality of the system is due to 
the physical origin of the cutoff. It appears as a consequence of the physical fact that far 
enough in the right spatial direction (the direction of propagation of the front) there must be 
no particles. This is, of course, totally independent of the spatial dimension of the system. 



10 



However, the effect of the cutoff on the dynamics of the front does strongly depends on 
the dimensionahty. In one dimension, we have observed a strong shift on the velocity of the 
front, while in two dimensions the effect is even stronger and the presence of the cutoff is the 
only responsible for the formation of diffusive instabilities jl^ . It would be very interesting 
to analize the effects of the cutoff in dimensions above d = 2, to see if new phenomenology 
develops or contrary there is a return to the mean field. 



V. CONCLUSIONS 



In this work, we have derived from first principles a field theory of reaction-diffusion 
particles, and we have used it to study reaction-diffusion propagating fronts. This kind of 
fronts has been traditionally studied using deterministic react ion- diffusion equations, like the 
Fisher equation, that consider an infinite number of particles. We performed a perturbation 
expansion in the ratio between the annihilation and the birth rates, that separates the mean 
field regime from the real discrete process, and studied the first order corrections to the 
equation of motion. A cutoff in the reaction term, a mechanism that has already been 
heuristically proposed, appeared in a natural way within our formalism, and leaded us to 
the velocity corrections already found in numerical simulations. 

Our first-principles analitically derived theory also allowed us to understand the funda- 
mental reasons that lead to the velocity shift. It is the compact support property of the 
front, i.e., the fact that the field is identically zero far enough in the spatial axis what causes 
such a dramatic effect in a pulled front like ours. 

It is also interesting to compare this work to a former one in the same direction Q| 
which tries to derive a cutoff in the reaction term, albeit for a different system. The main 
difference between both works is, under our point of view, that while this article concerns 
a system in the continuum space, the other deals with a lattice. This difference becomes 
fundamental since the calculations were performed using Stratonovich stochastic calculus, 
valid in the lattice, but ill-posed in a continuum space 16[. 

Many questions are still to be answered. Different reaction shemes are to be explored, 
also, the opposite limit (the annihilation rate large compared to the birth rate) is only 
conjectured jul], but not analitically found. Of course, higher dimensionality of the front 
is a very interesting problem, where new phenomenology does appear, as shown by Kessler 
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and Levine We hope that this and former works will encourage the reader to attempt 
to solve these and different problems that appear in this subject. 
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